Magnetic response and quantum critical behavior in the doped two-leg extended Hubbard ladder 
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We have investigated quantum critical behavior in the doped two-leg extended Hubbard ladder, by using 
a weak-coupling bosonization method. In the ground state, the dominant fluctuation changes from the con- 
ventional d-wave-like superconducting (SCd) state into density-wave states, with increasing nearest-neighbor 
repulsions and/or decreasing doping rate. The competition between the SCd state and the charge-density-wave 
state coexisting with the p-density-wave state becomes noticeable on the critical point, at which the gap for mag- 
netic excitations vanishes. Based on the Majorana-fermion description of the effective theory, we calculate the 
temperature dependence of the magnetic response such as the spin susceptibility and the NMR relaxation rate, 
which exhibit unusual properties due to two kinds of spin excitation modes. On the quantum critical point, the 
spin susceptibility shows paramagnetic behavior with logarithmic corrections and the NMR relaxation rate also 
exhibits anomalous power-law behavior. We discuss the commensurability effect due to the umklapp scattering 
and relevance to the two-leg ladder compounds Sri4_a;Caa;Cu2404i. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.10.Pm, 71.30.+h 



I. INTRODUCTION 

The two-leg ladder compounds, Sri4_2.Caa.Cu2404i, 
which exhibit the superconducting state' under pressure for 
a; > 12, have been studied intensively since superconduc- 
tivity was confirmed in a system with ladder structure. ^•^•"'•^•^ 
This material, consisting of chain and ladder layers, already 
has holes which are doped on the ladder layer even for the 
parent material, a; = 0. Reflecting its ladder structure, there 
exists a large gap in the magnetic excitationsii^ It is known 
that Ca substitution yields an increase of the doping rate from 
6 « 0.07 (a; = 0) to (5 w 0.25 {x = 12),'' and the material be- 
comes favorable for the superconducting state. A recent NMR 
measurement on superconducting materials under high pres- 
sure reveals evidence for a large spin gap (w 200 K)|i2iii4i^ 
which is much higher than the superconducting transition tem- 
perature Tc- From the measurement of the ^■^Cu NMR relax- 
ation rate T^^ at ladder sites, Fujiwara et al. '^ suggested that 
there are two excitation modes in the normal state: The one at 
higher temperature is an activation-type mode due to the spin 
gap and the other at lower temperature is an anomalous para- 
magnetic mode. In addition to the superconducting state, in- 
tensive studies have been devoted to the charge-density-wave 
(CDW) state, which is found in the parent system {x = 0) 
from optical measurements. '?il:ll Further the global phase 
diagram for overall hole dopingi^iii shows that, with increas- 



ing the hole doping, the CDW state is suppressed and dis- 
appears at a; ~ 9, while the superconducting state emerges 
for x ~ 11 and under pressure. From Raman scattering 
measurements,'^ it has been suggested that collective modes 
of the CDW state exist even in the highly doped superconduct- 
ing material a: = 12. Thus, it is expected that, by varying x 
and temperature, the competition between the superconduct- 
ing (SC) state and the CDW state will become crucial in these 
compounds, and nontrivial critical behavior will emerge in the 
competing region. 

Theoretical approaches to doped two-leg ladder systems 
have been performed by using various kinds of methods. It 
has been established that the d-wave-like superconducting 



(SCd) state becomes the most dominant fluctuation in the 
ground state of the doped two-leg Hubbard ladder and t-J 
ladder systemsi'^ffl^'i^V^i^^i^^i^^i^^i^^i^'^i^" When the model is 
extended to include parameters of several intersite Coulomb 
repulsions, other fluctuations can overcome the SCd state. A 
global phase diagram obtained from a weak-coupling g-ology 
approach shows that the CDW state, s-wave-singlet state, 
and d-density-wave states can also become quasi-long-range- 
ordered states in parameter spacei^Si^Si^ The CDW state has 
been confirmed by studying numerically the Hubbard model 
with nearest-neighbor repulsions. ^^ A richer phase diagram 
is obtained for the system at half filling^ ''^^ and at quarter 
filling. ^^"* Despite the huge number of theoretical works on 
the ordered state, the critical behavior expected close to the 
boundary between different phases is still unknown. In ad- 
dition, the temperature dependence of magnetic quantities, 
such as the uniform spin susceptibility and the NMR relax- 
ation rate, has been examined mainly on undoped spin ladder 
systems and the case with finite doping is not yet clarified. 

In the present paper, we investigate electronic states both at 
zero temperature and at finite temperature for the two-leg ex- 
tended Hubbard model with finite doping, where special atten- 
tion is focused on states near the phase boundary between the 
SCd state and the density-wave state. Based on the Majorana- 
fermion description, which is used for the low-energy effec- 
tive theory of the spin degrees of freedom, we demonstrate the 
unconventional temperature dependence of both the spin sus- 
ceptibility and the NMR relaxation rate close to the quantum 
critical region and also on the critical point, including loga- 
rithmic corrections. 

This paper is organized as follows. In Sec. |II] we intro- 
duce a model Hamiltonian for doped ladder systems and de- 
rive an effective theory describing low-energy physics, us- 
ing weak-coupling g-ology supplemented by bosonization and 
refermionization. In Sec. |III| we investigate the ground- 
state phase diagram by using the renormalization-group (RG) 
method, and clarify that the system exhibits a quantum critical 
behavior on the boundary between the SCd state and the co- 
existing state of CDW and p-density wave. In Sec. II VI the 



temperature dependence of both the spin susceptibiHty and 
the NMR relaxation rate is examined to clarify their anoma- 
lous behavior in the proximity to the quantum critical point. 
Finally in Sec. |V| we give a summary and discuss the com- 
mensurability effect due to umklapp scattering and the rele- 
vance of the present calculation to the experimental results for 

Sri4_a;Caa;CU2404i . 



II. MODEL HAMILTONIAN AND FORMULATION 

We consider a model of a doped two-leg Hubbard ladder 
with on-site and intersite Coulomb repulsions. The Hamilto- 
nian is given by 



H — Hq + Hu 



(2.1) 



The first term Hq describes electron hopping along and be- 
tween the legs: 



^o = -i|iE(4, 



l,<y ^j+i,l,<r 



j,<^,l 



^.E(4i 



<y j,2,o- 



H.c.) 



H.c), (2.2) 



The Fermi velocity of the bonding band and that of the an- 
tibonding band are given by vpfi = 2t|| cos{tt5/2 — A) and 
vf,tt = 2i|| cos(7r(5/2 + A), respectively. For (5 = 0, the Fermi 
velocity takes the common value vpfl — vf,tt — vp, where 
VF = 2i|| [1 - (i^/2t||)2]^/^ for arbitrary t^ (< 2i||). In the 
following, the difference between vf,q and VF,Tr is not taken 
into account since we restrict ourselves to the small-doping 
case \5\ <^ 1. 

Let us define order parameters of possible states. The most 
favorable state in doped ladders is the SCd state whose order 
parameter is given by^-' 



Oscd — 



j 



(cja,tCj,2,j, ~ c-,-44Cj,2,t) . 



(2.7) 



Other possible ground states are density-wave states with 
different angular momentai^ii^ In this paper we consider 
s-density-wave and p-density-wave (PDW) states. The s- 
density-wave state is nothing but the conventional CDW state. 
The order parameters of the density-wave states are given by 

Oa= J2 ^^(^) ^-(^) ^-('^ ± <5), (2.8) 

kM,± 



where Cjj^a annihilates an electron of spin a {—^, i) on the 
jth rung and ^th leg with I — 1,2. The Hamiltonian Hint 
represents interactions between electrons: 



Hint = U'^ nj,;_| rij^i^i + ^1 X! "j'' "j+i'' 



(2.3) 



where U (> 0), V|| (> 0), and V± (> 0) are coupling 
constants for the on-site repulsion, the nearest-neighbor re- 
pulsion on respective chains, and the nearest-neighbor re- 
pulsion on a rung, respectively. The density operators are 

Considering the two-particle interactions Hmt as a 
weak perturbation, we first diagonalize the single-particle 
hopping part Hq. The diagonalization can be per- 
formed in terms of the the Fourier transform, Ca-{k) = 
(27v)-i/2^^.^e-^'^"J-*'=^'cjj^^ where fc = (fc||, /c^) with 
k± — 0,Tr, and the lattice spacing a is set equal to 1. Then Hq 
is rewritten as Hq = J^k a ^i^) ^ti^) CCT(fe), where 



£(fe) 



-2i|| cos /cii — t± cos k±. 



(2.4) 



Here we consider the hole doping S for t± < 2i|| cos^(7r(5/2) 
where both the bonding (fcj^ — 0) and antibonding {k± = tt) 
energy bands are partially filled. In this case, the Fermi points 
are located at fc|| — ifcf.o and ±fci?,r for the bonding and 
antibonding bands, respectively, where 



kF,() ^ '^{l - S) + X, kF,^ = ^{l - 5) - X, 



and the quantity A is given by 



A 



sm 



t^/ [2^11 cos ^j] 



(2.5) 



(2.6) 



where A is CDW or PDW, and Q = (2fcF,7i") with kp = 
i7r(l — S). The form factors are given by /cdw = 1 and 
/pDW = isinA:||. These order parameters in real space are 
given by 

OcDW = 2j2{4,l,a'^J,l,<y - 4:2,^^.2,^) C0s2kFRj, 



3,0- 



(2.9a) 



OpDW ^ Yl (4-l,l,^S,l,<T - c]-l,2 



jA,cr'^j+l,l,<7 + Co,2,CTS+l:2,cr ) C0s2kFRj, 



(2.9b) 



where Rj = ja. The order parameter of the PDW state corre- 
sponds to that of the spin-Peierls state in the limit of 6 -^ 0. 
One can also consider other density-wave states such as the 
d-density-wave state, the /-density-wave state, or another su- 
perconducting state with s-wave symmetry, which are known 
to span the finite parameter space of the ground-state phase di- 
agram in the extended Hubbard ladder'"'^'' However, for the 
case of U ^ V\\, V± > 0, these unconventional states do not 
become dominant. 



A. g-ology 

Following the standard weak-coupling approach (g-ology), 
we linearize the energy bands around the Fermi points. The 
linearized kinetic energy is given by 

Ho= Y^ VF{pk\\ - kF,kJcl^^{k)Cp,^{k), (2.10) 



where the index p — + (— ) denotes the right- 
(left-)moving electron. We introduce field opera- 

tors of the right- and left-going electrons defined 

by ^p,.,+ (a:) = i"'/' Efc,, e*'^"""cp,.(fc||,0) and 

V'p,^,_(x) = i"'/'Efc|, e'''""cp,.(fc||,^) where L is 
the length of the chains: L = Na. 

The interactions near the Fermi points are written as Hint — 
J dx Hint, where 



^i- = lEE'[5il<.,a^- 



P,<T,C2 ^P,0-,C4 ^- 



s.'^'.Ca 



P,o- Ci = ± 



+ all %,,xi '^-P,ff.C2 ^P,ff,C4 V'-p.<.,c3 

+ 5211 ^laX, '/^-P,^,C2 ^-P,-,C4 V'p.a.Ca 

+ 52I V'p.^.Ci '^-p,aX2 ^-p,ff,C4 '/'p,<T,C3]' (2.11) 

where a =t (i) for cr =j (t) and e = C1C3 and e ee (1(2- 
The primed summation over Q (i = 1, . . . ,4) is taken under 
the condition C1C2C3C4 = +1^ which comes from the momen- 
tum conservation in the transverse direction. Each g'^'^ has two 
different processes [e.g., for 5++ one has (Ci, C2, Csi C4) — 
(+, +, +, +) and (— , — , — , — )], however; these two different 
processes are given in the same form in the bosonized Hamil- 
tonian and contribute to physical quantities in the same man- 
ner, as will be shown later. Then the coupling constants g^f. 
and gll are written in terms of interactions of the Hamilto- 
nian (12. 1> as 

4 = /,14+m,,,l^|, (2.12a) 

gti^U + kV^+m,,,V\\, (2.12b) 

where l± = ±1, toi.+ — — 2 cos tt (5 cos 2 A, toi._ = 
— 2cos7r5, 7712,+ = +2, and m2- — +2cos2A. 



B. Bosonization 

Here we apply the Abelian bosonization method i^^i^^iS The 
field operators of the right- and left-moving electrons are then 
written as 



V'p,a,c(2;) 



V2vra 



exp [ipkF,k^x + iptpp^s.ci^)] > (2-13) 



where s — + (— ) for a = f (J,). The chiral bosons obey the 
commutation relations [(pp,s.(;{x),^p.s'X'i^')] ~ ip7TSgn{x — 
x')ds,s'SQX' ^nd[ip+^sX^V-,s',C'] = «7r (5s, s"5(;,(;'. The Klein 
factors T^CT.c are introduced in order to retain the correct anti- 
commutation relations of the field operators between different 
spin and band indices. To relate the bosonic field ip to the 
physical quantity, we introduce a new set of bosonic fields 



bp± and (l)^± by 

(/>p+(a;) = ^ E 'Pp^-^xi^)^ 

sX=± 

(/>p_(a;) = J Yl C^P.^,c(2;) 



s,C=± 



- J2 '^fp,sxi^)^ 

sX=± 
sX=± 



(2.14a) 
(2.14b) 
(2.14c) 
(2.14d) 



The phases (j)p± and ipcr± represent charge and spin fluctu- 
ations, respectively, and the suffix + (— ) refers to the even 
(odd) sector 

In terms of bosonic fields, we can rewrite the kinetic energy 
as Hq — J dx Hq, where 



^o = ^EE[(^^ 



{dxcl^l 



u—p^(T r— ± 



(2.15) 



We also introduce the field ipi,r and its dual field 0i,r defined 
by 4i,yr = 4>tr + 4>ur and 0^^ = (t>tr ~ ^vr- Thcsc fields satisfy 
the commutation relation [01^^ (a;), 6*1/',' (x')] = — i7r0(— x + 
x')6r,r', where Q{x) is the Heaviside step function. 

In terms of these bosonic fields, the order parameters are 
expressed as Oa = J dx Oa, where 



Oscdix) 



pX 



Cipp,tci^)^-Pdxi^) 



oc e'^f+ cos 9p^ cos 



le 



suit 



sin (t>rr+ sin ( 



(2.16a) 



O. 



CDW 



O 



PDW 



(x) 



(x) 



P:0-,C 

OC sin (/>p_|_ cos 0p_ sin 00-+ sin 9^^ 



sinW„_ cos( 



E 

P,o-,C 



cos 6l<,_, (2.16b) 



+ sin0p+ sin0p_ cos(/)cr+ cosher-- (2.16c) 

In the present paper we do not consider the Akp CDW or- 
der parameter The possibility of the Akp CDW state will be 
discussed later in Sec.lIIII 

To obtain the bosonized Hamiltonian, Eq. (I2.13> is sub- 
stituted for the interaction term Eq. (I2.11> . The phase field 
(j)p- appears in the form cos(20p_ + 4Ax) where A is given 
by Eq. ( I2.6l l. We can safely assume that t± is relevant for 
t± being not very small. In this case we can discard the 
cos(2(/)p_ + 4Ax) terms which become irrelevant. We also ne- 
glect the cos 20O-- cos 26*0.- term, because this cannot become 
relevant due to the scaling dimension of cos 2(/)cr- cos 26* cr- 
being equal to or larger than 2, i.e., this term is either marginal 



TABLE I: Position of phase locking and signs for the fixed-point coupling constants, which is essentially the same as Ref.Bll The * symbol 
indicates that a bosonic field is not locked. liS are integers. 



W<T + 



-) 



-) 



^c-.s+' 



CDW+PDW (7r/2)(Jo + l)+7r/i {n/2)Io + Tvh * (7r/2)Jo + 71/3 (+,0,+,0,-) 

SCd (7r/2)Jo + ttJi (7r/2)Jo + vr/a (7r/2)Jo + 71/3 * (-,-,0,-,0) 



or irrelevant. Then our Hamiltonian reduces to 






r=± 



.P=± 



VF 

n 



E 



E (^-^'r 



9 = ± 



ttvf 



(9.0+ ) (5,0- 



g^^s+ cos26'p_ cos2(j)cr+ 
+ g^=^s- cos26'p_ cos20cr- 
+ ff^— COS 20p_ cos 26'o— 
+ 5s+.s- COS20O.+ cos20£r- 

+ 9s+.j^ cos 200-+ cos 29cr- 



(2.17) 



The coupling constants for the bilinear terms of the density 
operators are given by 



9p+ 


-+2>,(32t+52+r-5^ii), 


(2.18a) 


9p- 


= +\y]<9t\l+9tl-9{\), 


(2.18b) 


9<T+ 


-~\yM;-9ti~9n)^ 


(2.18c) 


9a- 


=-\y,<9t\i-9ti-9{\i 


(2.18d) 



9—^s+ 


= -9i± , 


ff^^s- 


= -9-,l. 


ff^— 


= +52I+ 


9s+,s- 


= +9tl, 


9s+.— 


= +5iT- 



and the coupling constants for the nonlinear terms are given 
by 

(2.19a) 
(2.19b) 
11+ , (2.19c) 

(2.19d) 
(2.19e) 

Since there is no cosine potential for the phase 0p+ in the 
Hamiltonian (I2.17> . the phase 0p_|_ is not locked even in the 
low-energy limit. Then the CDW state and the PDW state can- 
not be distinguished. Actually, by the translation of the phase 
(/)p+ — + (t>p+ +7r/2, the order parameters Ocdw [Eq. ( l2.16bH 
and OpDW [Eq. ( l2.16cH are interchanged, while the Hamil- 
tonian (I2.17> is invariant. Thus if the CDW state becomes 
(quasi-)long-range-ordered, the PDW state also becomes so, 
and then these two states coexist: We call this coexisting state 
the CDWh-PDW state. 



From Eqs. dzTel and dTTTt , the CDWh-PDW state and SCd 
state are identified by the fixed points of 5s as summarized in 
Tabled] Both states have a gap in the total spin sector 0^+. 
The first reason, attributable to the difference between the 
CDWh-PDW state and the SCd state, is that of the locking po- 
sition for the Op- mode and the 0^+ mode, where the solution 
gl_ ^^ < leads to {Op-) = ((/)^+) = or 7r/2 for the SCd 
state, and the solution g^ 4_ > ^ results in (0p_) 7^ {4'a+) 
for the CDW+PDW state. The second reason is the relevance 
of 0cr_ or 6'cr- : The CDW+PDW state is obtained for the lock- 
ing of 6'o._ due to the relevant g^ -— and g* -— terms, while 
the SCd state is obtained for that of (/)„- due to the relevant 
5^=_ and5*„_ terms. 



C. Refermionization and effective theory for spin modes 

The coupling constants in Eq. ( I2.11> are not independent 
parameters due to the global spin-rotation SU(2) symmetry. 
In terms of the coupling constants in Eq. ilAli . the constraint 
is given bjJ^^i^ 

5(T+ + 9<j- - gs+,s- = 0, (2.20a) 

9a+ - 9cr- - 9s+.— = 0, (2.20b) 

9^.s+ - 5— s- - 5—— = 0. (2.20c) 

Since the SU(2) symmetry holds in the original Hubbard 
Hamiltonian ( 12. 1> . the coupling constants in Eq. (I2.17> must 
satisfy Eq. ( I2.20> in the course of renormalization. 

To appreciate the SU(2) symmetry in the effective theory 
(I2.17> . we fermionize it by introducing spinless fermion fields 
'0p,r (p = ± and r — ±): 



V'±,r(a;) 



Vr 



exp[±i20+^(x)] 



(2.21) 



where the index r = +a (— ) refers to the total (relative) de- 
grees of freedom of the spin mode, and {rir,rir>} — 2Sr,r'- 
The density operators are given by : ■01 ± V'p ± • = dx(j)^±/Tr. 
We then introduce the Majorana fermions ^" (n = 1 — 4) by 



V'p,+ = 7f (^P + *^p)' ^P'- " 7f (^: 



^O- (2.22) 



These fields satisfy the anticommutation relations 

Up(^),Q'(^')} = S{x ^_^Spydn,n'- With the help 



of the SU(2) constraints ( I2.20> . we rewrite the effective 



Hamiltonian (I2.17> in terms of the Majorana fermions: 



n 



Vp 



jp Y^ 



E 



'pr) 



9pr 

TTVp 



(a,0+ ) (a,0^ 



pr/ 



■ ^C — ,St (-1 /I 



i^+-i-y 



27ra 

. yc—,ss 
i—— — 

27ra 



cos26ip_ C+-C- 



-9.-{i+-i-)c+e-, 



(2.23) 



where ^^ = (^p, Cpj Cp)' ^"d the couphng constants are 

9—M = ^ff^,s+^ 9—,ss = -9^,s- + ff^,—- (2-24) 

These coupUng constants are given in terms of the Hubbard 
interactions as 

5— .si ^+U-Vj_~ 2V\\ cos ttS, (2.25a) 

5— ss ^+U -V^ + 2V||(cos7r5 + 2cos2A), (2.25b) 
5p+ = +U + 2Kl + ^i [4 + cos7r(5(l + cos 2 A)], (2.25c) 
5p_ = -Fl - V|i cos7r(5(l - cos2A), (2.25d) 

5„+ = +;/- yj|COS7ri5(l + cos2A), (2.25e) 

5„_ =+Kl + V||COS7r(5(l-cos2A). (2.25f) 

Thus the effective theory for the spin sector becomes 0(3) x Z2 
symmetric, i.e., the four Majorana fermions are grouped into a 
singlet ^'' and a triplet ^. We note that the 0(3)xZ2 symmetry 
also appears in the low-energy effective theory of the isotropic 
Heisenberg ladderiS 

III. PHASE DIAGRAM IN THE GROUND STATE 

We investigate the low-energy behavior using perturbative 
RG analysis. There are six independent RG equations for the 
scaling of the coupling constants under the transformation of 

From Eq. J2.23> . we obtain the 



the lattice constant a 


-^ ae^'. 


RG equations 




s«'--i-fc 


,st ~ ^^ 


-G,+ ^ -Gl^ - 


-Gl_- 


-G,_ = -2G,+ 


G,_- 


d^ r 


^■'€-,31 


d^ r 


- G'c-.ss 



y2_ 
c—,ss ' 



2 c-,st' 

_/^ (^ 

r\ C—.St C — SS^ 

- 2G0-+ G— f 



oG(j— G~ 





(3.1a) 




(3.1b) 




(3.1c) 


G,_ 


-Gc-,ss^ 




(3. Id) 




(3.1e) 



and dGp+/dl = where G(0) = g/{2Trvp). Note that these 
RG equations can also be derived directly from Eqs. (I2.17> 



and J2.20> . Since the coupling Gp+ is unchanged under renor- 
malization, the total charge sector is critical and has gapless 
excitations. The asymptotic behavior of these coupling con- 
stants for large I is examined by integrating the RG equations 
(13. 1> numerically with the initial conditions ( I2.25> . It is eas- 
ily found that, in most cases, all the coupling constants in 
Eq. ( 13. 1> grow under renormalization and become relevant at 
large I. This fact implies that the all the modes except for 
the (j)p+ mode become massive in most regions of the ground- 
state phase diagrami^2i22iiii These stable fixed points are called 
the "CI SO phase," where the notation CnSm denotes n mass- 
less boson modes in the charge sector and m massless bo- 
son modes in the spin sector ^^ The characteristic energy scale 
corresponding the mass gap can be roughly estimated from 
I nia I — A e^'" where A is the high-energy cutoff of the order 
of the bandwidth and la is determined by using the fact that 
the corresponding coupling constant |Ga(OI becomes of the 
order of unity at Z = /„. 

From Eq. ( 13. 1> . we also find that there are two distinct sta- 
ble fixed points on the plane of the coupling constants for the 
original extended Hubbard ladder model with U (> 0), V|| 
(> 0), and Vj_ (> 0). Actually, from Table 12 one obtains 
that the coupling constants {Gp- , Gc+ ,Ga-, G— ^^ , G— ^j, ) 
flow to (-,-,+,-,+) for the CDWh-PDW' state and 
(— , — , — , +, +) for the SCd state. This means that there are 
two distinct phases in the ground states of the extended Hub- 
bard ladder model, and that the system exhibits a quantum 
phase transition on a critical point between two phases. In or- 
der to examine critical properties in the ground state, we ana- 
lyze Eq. ( I2.23> in more detail by deriving the effective theory 
for low-energy properties. 

Here we assume that the mass of the charge mode of odd 
sector (p— ) is larger than those of the spin modes {a±), so that 
the 9p- fields are locked by cosine potential below the scale 
of the mass rup- . This assumption will be examined later 
The effective low-energy theory is obtained from Eq. 
by taking an average: 



(cos26'p_). 



(3.2) 



Then we have 

Vp 



|_-a,^_)-zm?^+.|_ 



*f (eta.a-e5.e)-^m°ete 



y<T+ 



{£.+ ■£.-) -9.- {i+-i-)e+ it (33) 



where we have introduced 



m° 



:ffc 



:3c 



(3.4) 



Zna ' lira 

Such a mean-field treatment of the charge sector has also been 
utilized in the context of carbon nanotubesill We note that 
this low-energy effective theory takes the same form as that 
of the isotropic Heisenberg ladder.^^ In terms of the original 
Hubbard interactions these masses of Eq. \2>A\ are given by 



m° = -^ [[/ - Vl - 2V|| cos7r5] , 



(3.5a) 



m" 



p- 
2'na 



[[/ - VI + 2V|| (cos 77^ + 2 cos 2A)] . (3.5b) 



Thus the magnitude of the masses can be tuned by the inter- 
actions and doping. It is known that, when m", m'^ ^ 0, the 
quartic marginal terms lead to mass renormalization, m° -^ 
nis and to? --> rrit, whera^i^ 







mt — rrif 



TO, 



-TOj In -—Q 



TTVp 



2t:vp 



m^ In 



A 



r, (3.6a) 



(3.6b) 



2ttvf i'"-*! 

and A is a high-energy cutoff. Then Eq. J3.3t reduces to 

Low-energy properties become more transparent by intro- 
ducing four copies of the one-dimensional quantum Ising 
model: 



HQi = -l^l^{J-h-. 



]+id 



hl<l) 



(3.8) 



where a^ and a^ are the Pauli matrices and / = 1,2,3,4. This 
model is equivalent to the Majorana-fermion theory with the 
central charge c = 1/2. The operator pi, being dual to the spin 
operator a, is defined as^ 



M 



'j+l/2,l 



H' 



;,i' Mj+i/2,i - ^j,i^j+i,i> 



(3.9) 



which is known as the Kramers-Wannier transformation in the 
one-dimensional quantum Ising model. These variables sat- 
isfy [<,,, Mj+i/2,;] = for i > j and {<,, M^+i/a,;} = for 
i ^ j- In terms of a^ and /i^, i.e., the Ising order and disorder 
parameters, the Majorana fermions can be constructed as 

where k; is the Klein factor. One can easily check 
the anticommutation relation of the Majorana fermions, 
{Vi.hVj,™.} = {Ci,hC]..m} = 2SijSi^m and {■r]i,iXj,m} = 0. 
By using ^i = {-m + Ci)/V2 and ^i = {rji + Ci)/V2, and 
by taking the continuum limit, the quantum Ising Hamilto- 
nian Eq. ( 13. 8> reproduces Eq. (13. 7> where vp — 2 J, nit — 
2{hi — J), and to^ = 2(/i4 — J) with hi = h2 = hj,. It 
is well known that the Ising model ( 13.81 1 exhibits a quantum 
critical point at hi = J.*^ For hi < J, the ordered state 
is obtained, i.e., the order parameter ai has a finite expec- 
tation value. For hi > J, on the other hand, we have the 
disordered state where the expectation value of cri becomes 
zero, while the disorder parameter /i; has a finite expectation 
value. On the critical point hi = J, the corresponding mass 
in Eq. ( 13.71 1 vanishes with its central charge c = i for each 
Ising chain. '^^ Thus the ground-state properties are determined 
from the sign of masses toj and m^. When mt < 0, i.e., the 
Ising model with I = 1, 2, 3 is in the ordered phase, we have 
(cti> = (0-2) = (0-3) ¥" and (/ii) = (^2) = (ms) = 0, 



TABLE II: Possible phases and related quantities: the signs of 
masses (mt and ms) and order parameters. We have assumed 



c—={6,.)^ 


: mod TT in Eq. <3.4t. 


rrit nis Order parameters 


CDW -1- PDW 

SCd 


+ (^1,2,3) / 0, (^4) / 

+ + (M1,2,3) / 0, (m4> / 



and vice versa. In the same manner, we have (0-4 ) 7^ and 
{114) = for nis < 0, while ((74) — and {114) ^ for 
TOs > 0. In terms of the Ising variables, the order parameters 
Eq. (I2.16t are rewritten as 

Cscrf«e*^''+ (/^i^2M3)Ai4, (3.11a) 

CcDW OC Sin0p+ (criCr2CT3)A^4, (3.11b) 

OpDW OC COS 0p+ (criCT2 0-3)^^4 • (3.1 Ic) 

As noted in the preceding section, we find that both the CDW 
and PDW states have the same structure for the spin degrees 
of freedom, since the field (\>pj^ is unlocked due to the doping 
effect. When to* < and m^ > 0, the CDWh-PDW state 
becomes quasi-long-range ordered. In the case m^ > and 
TOs > 0, the dominant fluctuation is the SCd state, which 
is called the "Luther-Emery liquid.'"*^'*^ The possible ground 
states and those order parameters are summarized in Table HIl 
Let us examine the behavior in more detail using the scaling 
equations for the coupling constants in the effective Hamilto- 
nian ( 13. 3> . The scaling equations for the coupling constants 
are given byS 



dMt 
dl 
dMs 
"dF 
dG^ 

dl 
dG^ 

dl 



^Mt-2MtG^+-M,Ga-, (3.12a) 

= Ms-3MtG,-, (3.12b) 

- = -Gl+-Gl_-M^, (3.12c) 

- = -2G,+G,- - MtM„ (3.12d) 



where dl = da/a, Alt — m^a/vp, Ms — m'^a/vp, and 
Ga-± = ga±/2Trvp. The couplings Mt and Mg are relevant, 
while Gcr± are marginal. These RG equations, which are an- 
alyzed in a way similar to Ref. 33, have two kinds of stable 
fixed point M^ = ±cx3. The ground-state phase diagram, 
which is summarized in Table HIl is shown in Fig. Q] on the 
plane of U/t and V/t (ty = tj. = i, V|| = Kl = V). The SCd 
quasi-long-range-ordered state, which is obtained for V — 0, 
is destabilized by the intersite Coulomb repulsion and changes 
into the CDWh-PDW state. The ground-state phase diagram 
on the plane of the doping 5 and the ratio V/U is shown in 
Fig. 12 where the SCd state is stabilized due to the suppres- 
sion of the intersite repulsion by the doping. 

Here we discuss the possibility of the Akp CDW state. It 
is known that the Akp CDW state also becomes quasi-long- 
range ordered in the whole region of the phase diagrams. Fig- 
ures ^ and |2] since the correlation function Akp CDW state 
decays^^'^*^ as l/r'^^p+ where Kp+ is the Tomonaga-Luttinger 




u/t 



FIG. 1: The ground-state phase diagram on the plane of U/t and 
V/t where t = t± = iy and S = 0.2. 
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parameter for the total charge sector, Kp+ = [{2Trvp — 
gp+)/{2'KVF + 9p+)V^^- Since the correlation function of 
the CDW and PDW states decays as l/r^p+/^, we find that 
the Akp CDW state is still a subdominant fluctuation in the 
CDW+PDW state. On the other hand, in the SCd state, the 
ikp CDW can become the dominant fluctuation for Kp^ < 
1/2, since the exponent of the SCd correlation function is 
given by l/2Kp^. However, at close to half filling, it has been 
confirmed that the exponent Kp+ in ladder systems reaches 
universal value K*^ -> 1 as (5 ^ 04^'*^''*** Thus in the SCd 
state of the phase diagram, we find that the correlation func- 
tion of the SCd state becomes dominant and that of the Akp 
CDW state is subdominant close to half filling. 

The doping dependences of the gapsp^ which are roughly 
estimated from \ma\ — Ae~'", are shown in Fig. |3l The 
Majorana triplet gap nit collapses on the boundary between 
the CDWh-PDW state and the SCd state where the sys- 
tem is strongly fluctuating due to the competition between 
these two states. From the perturbative RG method, we 
cannot determine the precise magnitude of the gaps; how- 




0.4 0.5 

doping 5 

FIG. 2: The ground-state phase diagram on the plane of doping rate 5 
and ratio V/U where V = V« ^V± and U/t = 3 with t = t±^tu. 



FIG. 3: The doping dependence of the energy gaps nit and m^ with 
U/t = 3, V||/f = V±/t = 0.7, and i = ty = tx = 1. The quantum 
critical point (QCP) is at <5 = 5c « 0.18. The CDW+PDW quasi- 
long-range-ordered state is obtained for 5 < Sc, while the the d-wave 
SC quasi-long-range-ordered state is obtained for 5 > 5c. 



ever, the qualitative features of the gap associated with the 
phase transition (i.e., mt in the present case) can be captured 
close to the quantum critical point (QCP). The critical prop- 
erties are described by the SU(2)2 Wess-Zumino-Novikov- 
Witten (WZNW) model"-^^^^"^. and is also characterized by 
the C1S| stateii^ In the present analysis, we have assumed 
TOp- > {nis, \mt\) to derive the effective low-energy theory 
Eq. (13. 3> . Actually, by using Eq. ( 13. 1> . we find that nip^ be- 
comes the largest among the three nip-, mt, and to^. How- 
ever, nip- is not much larger than the Majorana singlet gap 
nis, but is of the same order as rn^; we find from the above 
rough estimation that r7ip_ /tos w 1.4— 1.6 for 5 = 0.1 — 0.4. 
From the recent theoretical studies on the multi component 
one-dimensional systems, it has been proposed that^^-^-'-^'* a 
symmetry, which is broken even in the microscopic Hamil- 
tonian, can be restored nontrivially at low energies, due to 
coupling terms between different modes. This mechanism 
is called a dynamical symmetry enlargement (DSE) whose 
possibility has also been examined by the nonperturbative 
approachi^ Once the DSE is realized, the resultant gaps for 
different modes become identical. Based on these theories. 



the present results with nip 



may suggest the occur- 



rence of the DSE between the charge sector (p— ) and the Ma- 
jorana singlet sector where excitations form an 0(3) multiplet. 
However, the study of the DSE is beyond the naive pertur- 
bative RG approach^^ and thus the DSE in the present case 
remains unclear. 



IV. SPIN SUSCEPTIBILITY AND NMR RELAXATION 
RATE 

In this section we study the uniform spin susceptibility and 
the NMR relaxation rate at finite temperature from two differ- 
ent approaches by extending previous calculations on the sin- 
gle chain''^-'''' or on the undoped Heisenberg 1 adder ^IsSi^Siffi 
One approach is the random-phase approximation (RPA) 
combined with the renormalization group method and an- 
other is direct calculation in terms of the low-energy effective 
Hamiltonian (13. 7> . The former has the advantage of reproduc- 
ing high-temperature behavior, while the latter is appropriate 
for describing the low-temperature asymptotics. 

Fkst we introduce the spin-i operator; 

^(9) = ^ E 4(fc + '7)T,,,..c,.(fc), (4.1) 

fc,(Ti,CT2 

with q — ((7||, (7^) and the Pauli matrices (Tcr^^a2- Ths gener- 
alized spin susceptibility is given by 

X{q, liUn) ^^J dT {TrS'^iq, r) 5"(-g, 0)) e"^"\ 

(4.2) 
Equation ( 14. 2> is independent of a due to the spin-rotational 
SU(2) symmetry, where a stands for the orientation of the 
magnetic field. The noninteracting susceptibility [i.e., Eq. 
( 14. 2> without interactions] is given by 



Xoiq,iuJn) 






f[e{k + q))-!{e{k)) 



4L ,^ iujn + e{k)-e{k + q)' 



(4.3) 



where /(s) is the Fermi distribution function /(e) = 
l/[e'^'^'^~''-' + 1], and a;„ and /i are the Matsubara frequency 
and the chemical potential, respectively. In the continuum 
limit, we can split the spin operator into a uniform part varying 
slowly in space and a staggered oscillation part, as 



S(x,q^) ^ Jr{x) + i-lf/'^nrix), 



(4.4) 



where r = + (— ) for gj_ = (vr). The uniform part {q\\ « 0) 
of the spin operator is given by 






<Ti,^2V'p,a2,c(^)' (4.5a) 



P,C f7-l,CT2 



where V'p.o-.cl^) is given by Eq. (I2.13> . By using Eqs. ( I2.13> . 
( I2.21> . and ( I2.22> . the spin operator J±(x) is expressed in 
terms of the Majorana fermions as 



J+{x) = +-^ip{x)y.$,j,{x), (4.6a) 

p 

J-{x)=-iY,ip{x)ep{x). (4.6b) 
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FIG. 4: The spin susceptibility x(9||,?±,^) with gx = (upper 
diagram) and tt (lower diagram), and small q\\ . The subscript p = + 
(— ) refers to right- (left-)moving electrons, while C = + ( — ) refers 
to electrons on the bonding (antibonding) band. 

where J± = {J^, J^, J^) — (J4,J±,J±). The staggered 
part ((7|| « 2kF and q± = tt) of the spin operator is given by 

"-(2^) = 2 ^ ^ ^Li,c(^)'^-i.-2V'-P,-2.-c(^)- 

P,C ''■l:'''2 

(4.7) 
Here we do not consider the component with {q^^,q±) w 
{2kp,Q), which would become irrelevant in the low-energy 
limit due to the relevant t±.^^ By using Eq. ( I2.13t . the 
tor n. is rewritten as 



t in the low-energy 
q. ( I2.13t . the opera- 



-2i 



— - ZZ " 

nt ix) — (cos 0P+ cos 0p_ sin 6*^+ cos (})„- 

na 



— sin </)p_|_ sin Op- cos 0^+ sin (f>cr-), (4.8a) 



-2i 



—2i 

Tv!_ (x) = (cos 0p+ cos 9p- cos 9a-+ cos (j)c 

na 



cos (j)a- 



sin (/)p_|_ sin 6'p_ sin 9a-+ sin (pa--), (4.8b) 



2 

nl (x) = — (cos 0p+ cos 0p_ cos (/)„+ sin 0cr- 
7ra 

— sin (f>pjf- sin 0p_ sin 4)^+ cos 9a-)- 

c. Here we note that n_ (a;) expresses 



— sm cppjf- sm Wp- sm (^ 

where 0p+ = (/)p+ — ttSx. Here we note 
the incommensurate spin-density wave. 



(4.8c) 



A. Uniform spin susceptibility 

Here we calculate the uniform spin susceptibility. In terms 
of the spin operators J±, the spin susceptibility x{q, i^n) for 
small q|| is given by 

1 f^ 
Xuni(g,*c^„) = ^j dr(T,J,"((j||,T) J,"(-g||,0))e''^"^ 

(4.9) 
where r = + (— ) in the right-hand side (RHS) corresponds 
to g^ = (tt) in the LHS and a denotes the orientation of the 
magnetic field. First we calculate Xuni(Q, i^n) within the RPA 
by using the formulation of 5-ology ( 12. lU . The magnetic sus- 
ceptibility x{% 5 1^1 i^n) for small qn is calculated from the 



diagrams of the RPA given in Fig. |4] We note that the same 
result would be derived from the path integral formalismJ^ 



The explicit forms for the susceptibility are given by 



Xuni{q\\,0,tujn)^ ++ __ + ^ -— , (4.10a) 

1- [giZ +9i± ) Xo('?||,0,it^„)xo [Q\\,0,tuj„) 

I 



where Xo(9|| j 1^ ■• *^") i^ given by the noninteracting spin sus- 
ceptibility per branch p = + (— ): 



Xg(q,i^„) = — Y^ Y^ 



f{e{k + q))-f{e{k)) 
iujn + s{k) — e{k ~\- q) ' 



ku^pkjp k_[ 

(4.11) 
It is crucial to take into account the effect of the curvature 
of the dispersion in the noninteracting susceptibility Xo' i-^-, 
the explicit form of e{k) given in Eq. (12. 4> is retained in the 
calculation, to obtain a reasonable temperature dependence. 

By using Eq. (I4.10a> with the relations Eqs. (I2.19> 
and J2.20L the uniform spin susceptibility Xs(T) — 
Xuni(g||, 0,0)1^11^0 is calculated as 



Xs{T) 



Xo{T) 



1 -.9<T+xo(r)' 



(4.12) 



where xct(0,0) = Xcr(0,0) [=xo(T)/2]. We note that the 
susceptibility Xo{T) for T = is given by Xo(0) = {vp.o + 
VF.7r)/{4:T!'VF,QVF.Tr), which rcduccs to xo(0) -^ 1/{2ttvf) in 
the limit of a single chain because vf,q, vf,-r — > vf- 

All the diagrams given in Fig. 0] are non singular because 
only the small-^n components are taken into account. The ef- 
fect of the one-dimensional fluctuations appears through loga- 
rithmic corrections to §,,+ . This process can be accomplished 
by replacing the coupling constant 9^+ with the renormal- 
ized one with the cutoff for finite temperature, i.e., ga+{l — 
In A/r). Then the uniform spin susceptibility is given by 



Xs{T) 



Xo{T) 



i-5.+(0xo(r)' 



(4.13) 



where ga+{l) [= 27rwi?Gcr+(0] is obtained by solving the RG 
equation Eq. ( I3.1b> or ( I3.12c> . This formula is valid at the 
temperature above the energy scale of the gap \int\, since the 
coupling (7ct+ which describe the fluctuation of the Majorana 
triplet sector [see Eq. (13. 3H is treated perturbatively. 

The overall temperature dependence of the uniform spin 
susceptibility is shown in Fig.|5] At sufficiently high temper- 
ature, the uniform susceptibility Xs (T) exhibits behavior sim- 
ilar to the one-dimensional susceptibility both in the presence 
and in the absence of interactions. With decreasing temper- 
ature, Xs{T) and Xo{T) for the ladder system become larger 
than those of a single chain due to the enhancement of the den- 
sity of states by interchain hopping. However, Xs{T) shows 



activation behavior below \int\. In order to comprehend such 
low-temperature behavior, we estimate the spin susceptibil- 
ity based on the effective theory for the spin mode, i.e., the 
Majorana-fermion theory. Below the energy scale of the gap 
\mt\, we could ignore the fluctuation effects due to go-+ in 
Eq. ( 13. 3> . which would merely yield the mass renormaliza- 
tion given in Eq. (I3.6a> . Thus we can use the effective theory 
given by Eq. ( 13. 7> . although there remain some discussions on 
ga+-^'^ By using Eq. ( 14. 6> and after a straightforward calcula- 
tion, we obtain the uniform susceptibility as 



Xs{T) = XuniCgil, 0,0)1^1,^0 



1 



dk sech 



where e J, = (vpk'^ - 



'k - y"F 

is rewritten as 



STrry 2T' 

(4.14) 
1/2. For low T (< \mt\), Eq. (l4TT4b 



1 



-m^^VS-'-"^' 



(4.15) 



Thus we obtain the exponential decay of the spin suscep- 
tibility in the doped Hubbard ladder, which is the same 
low-temperature asymptotics as in the undoped Heisenberg 
ladder JiS 

The susceptibility at the QCP exhibits quite different be- 
havior from that of y = 0. Figure |6] shows the ladder Xs{T) 



Q 

IS 

•S 0.08 




' '^^^^^'^^''^^^^ 




^:::= 


=^ 


S 




ladder Xs 


- 


tZ3 

0.04 




ladder %» 

single chain Xs 

single chain X(i 


- 





, 







Temperature T 

FIG. 5: The temperature dependence of several spin susceptibilities 
Xs{T) and xo{T) for U/t = 4, V^ = V± = 0, and 5 = 0.1, with 
i|| — t± — 1, where xo denotes Xs in the absence of interactions. 
The Majorana triplet gap is nit ~ 0.03. For comparison, the corre- 
sponding susceptibilities of a single chain are also shown. 
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FIG. 6: The temperature dependence of the spin susceptibility for 

U/t = 3, V|| = VI = 0.7 with f y = tx = 1. The critical value of 
the doping is given by 5c ~ 0.18. 



at low temperature for the doping near the QCP. Except for 
the QCP, the susceptibility exhibits activation behavior at low 
temperature. On the QCP, the susceptibility shows paramag- 
netic temperature dependence for the whole temperature re- 
gion due to the absence of the gap rrit- 



B. NMR relaxation rate 

In this section, we calculate the temperature dependence of 
the NMR relaxation rate T^^. We use the general formula of 
the NMR relaxation rata 



,M 



1 2L^ uuo ' 



(4.16) 



where x" is the imaginary part of the magnetic susceptibil- 
ity and ijjQ is the nuclear resonance frequency having a small 
energy of the order of millikelvins. Here we do not write the 
coefficient of the RHS of Eq. ( I4.16> by neglecting the momen- 
tum dependence of the hyperfine coupling constant, since the 
hyperfine coupling for the ladder-^^Cu site in Sri4Cu2404i 
originates mainly in the on-site hyperfine interaction and thus 
depends weakly on the momentum q? Since the interaction 
process between electrons near Fermi points is considered to 
play an important role for weak coupling, the integral over q\\ 
in Eq. ( I4.16> can safely be split into two parts. 



T^'^{T^'U, + {T^') 



stag! 



(4.17) 



where {T-^ ^)uni and [T-^ ^)stag are the contributions from 
small g|| and large q\\, respectively. Thus (Tj"^)uni denotes 
the uniform contribution to the NMR relaxation rate, while 
^stag is the antiferromagnetic staggered contribution. 



The explicit formulas are given by 



{T^\ 



T 



-y 



An 



dq\ 



Xu„i(9'^o) 



q± •"'II 



911 ~0 



Wo 






x"iQ,^o) 



dq\ 



(4.18a) 



(4.18b) 



These two contributions (Tj~^)u,ii and (Tj~ )stag show dif- 
ferent temperature dependence. It has been discussed 
which contribution becomes dominant in Heisenberg ladder 



systems? 



57.58.64.65 



but the problem is still controversial. We 



treat these two contributions separately in the following. 

1. Uniform part: Contribution from small qi^\ 

Now we examine the uniform contribution (rf ^)uni, by us- 
ing the RPA calculation combined with the RG method at high 
temperature and by performing a direct calculation in terms of 
the low-energy effective theory at low temperature. 

First we focus on the high-temperature behavior of 
(Tj~^)uni- In order to obtain (Tj~^)uni from Eq. ( I4.18a> . we 
calculate the imaginary part of the susceptibility Xuni- ^^ the 
noninteracting case, one easily finds from Eq. ( 14.1 1> that the 
imaginary part of the noninteracting susceptibility becomes 

Xo('?|hO,w) = -^pq\\5{uj-pvFq\\) (4.19a) 

p 

Xo(9|h7r,Lj) = -'^{pq\\+2C\)6{uj-pvpq\\~2CvF\), 

(4.19b) 



p,C 



where Xo (g|| , g_L, t^) = Ep=± 1^X0(911 , ^.l, ^) and we have 
neglected the curvature of the dispersion and used the Un- 
earized dispersion e(fc||,fc^) -^ vp{pk\\ — kp.k^) in Eq. 
(I4.11> . Since Xo (9i ^) is proportional to u) due to the 5 func- 
tion, the integral in Eq. ( I4.18al l is determined by the contri- 
bution being linear in Xo in Xuni((Z|| , 0, uf) [Eq. ( l4.10aH . Fur- 
ther, by neglecting the q\\ and u dependence of Rexo(9|| : 0, oj) 
(which is a nonsingular quantity) in Eq. j4.10aL the imaginary 
part of Xuni(9|| , 0, lo) for small q\\ and small lo is given by 



Xu„i(9||,0,w) 



Xo(9|hO,tj) 



[l-5,+x(,(0,0,0)]' 



4x§(r) 



^ pq\\ 5{uj - pvpqw), (4.20) 



where Xo(9||,9_L,t^) = Ep^± ReXo(g|| . g-L, t ^) and xo(T) ee 
x'(0,0, 0). In the second equality of Eq. ( I4.20> . we have 
used Eqs. (I4.13> and ( I4.19ak On the other hand, there is no 
correction to the q^_ — n component, i.e., Xuni(g|h "'' ^) ~ 
Xo(<7||,7r, w). By inserting Eqs. ( I4.20> and ( I4.19b> into Eq. 
(I4.18a> . the uniform contribution (T]"^)uni at T ^ m^, |mf | 
is given by 



(Jl juni — 



T xM. 
xl{T) 



T 



SttvJ, 



SttvJ,' 



(4.21) 



The first (second) term in the RHS of Eq. ( 14.2 1> comes from 
processes with momentum transfer q± = (q± = tt). The 
contribution with q± — tt, which takes the same form as in the 
noninteracting case, shows the Korringalaw T-^^ ex T, while 
the contribution with g^ = is enhanced by the spin fluctua- 
tions yielding a factor xliT)/xl = 1/[1 - 9a+{T)xo{T)]^- 
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The contribution of q± = shows a relation between the uni- 
form part of the relaxation rate and the uniform spin suscepti- 
bilityx.(r). 



Next we calculate (1/Ti)u,ii at low temperature by using 
the low-energy effective spin Hamiltonian J3.7> which is valid 
at T <C TOp_. By using Eqs. (13. 7> and (14. 6> . we obtain 



J 



{Tr'). 



Wnvp 



de —j= 

OO 



n[(£ 



-Wo J 

-2 I 



de ■ 

V 



■mtrris 



sech 



167ri)^ 7max(|mt|.me) ^J {s^ - mDiiz^ - mf) 



2T 



sech^— -, 



2T' 



(4.22) 



{vj 



ly^'^ with z — t,s and we have set 



Wo — > in the second term of the RHS. The first term in the 
RHS of Eq. ( 14.221 1 is a contribution from the processes with 
momentum transfer q± = 0, which is given by the Majorana 
triplet-triplet bubble in the diagram. ^^ On the other hand, the 
second term is a contribution from processes with q± = tt 
and corresponds to the Majorana triplet-singlet bubble. In the 
limit of high temperature, Eq. ( I4.22t reduces to (Tj^^)uni — 
T/Snvp + T/SiTVp corresponding to the Korringa law. The 
enhancement factor in the first term of the RHS of Eq. (14.2 1> 
is not reproduced since the spin fluctuation has been neglected 
in Eq. ( I4.22> . At low temperature (wo -^T -^ |"it |), the most 
dominant term reads 



{Tr'y 



Anv'j, 



de- 



(£2 + m2)e-^/^ 






invp 



^t\ ^-Wt\IT 



AlTVp 



\2TJ 
Infill -7 



(4.23) 
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FIG. 7: The temperature dependence of the uniform contribution 
of the NMR relaxation rate, with U jt = 3, V\\lt = Vx/t = 0.7, 
t = t|l = t^ = 1, and 5 = 0.1. The dashed and dotted curves 
denote the gx = and -k contributions, respectively, and the solid 
curve denotes the total (Tj^^)uni. 



where Kq[z) is the modified Bessel function and 7 is Euler's 
constant. Equation ( I4.23> is equal to the formula obtained in 
the undoped two-leg Heisenberg ladderiSi^ and in the spin- 1 
Haldane spin chain. ^i 

The overall temperature dependence of (Tj~^)uni can be ob- 
tained by the interpolation between Eqs. (I4.21> and ( I4.22> . 
We show the temperature dependence of (Tj"^)uni in the 
CDWh-PDW phase in Fig. In the SCd phase, (Tf ^)uni 
shows behavior qualitatively similar to Fig. The dashed 
curve represents the gx = contribution while the dotted 
curve represents the q^=-K contribution. The q_L = tt contri- 
bution exhibits activation behavior at temperatures below m^, 
while the gj_ = contribution shows T-linear dependence at 
temperatures above \mt\. For T <C |mf |, the NMR relaxation 
rate (rj~^)uni is governed by the first term in Eq. ( I4.22> and 
exhibits an activation behavior given by Eq. ( I4.23> . The com- 
ponent with qx = is always larger than that with q^= ixm 
the whole temperature region, in contrast to the results in Ref. 
dTI It is found that the total (Tj^^)uni clearly exhibits behav- 
ior consisting of two components and thus the result cannot 
be fitted by a single activation energy. 



2. Staggered part: Contribution from large q^^ 

We calculate the antiferromagnetic contribution {T^^)sta.g 
given by Eq. ( I4.18b> . As is known in the single-chain casei^ 
both the spin and charge degrees of freedom contribute to 
(Tj^^)stag- The critical charge mode in our theory is the total 



sector 



, whose Hamiltonian is given by 



ZTT 



K. 



p+ 



-[d,<i)p+Y + Kp+{d^ep+Y 



, (4.24) 



where Vp+ = vp[l — {gpj^/2'KVpYY/'^ and Kpj^ — [{2'iTVp — 
gp+)/{2'KVF + gp+)Y^^ ~ 1 — Gp+. For simphcity, we will 
neglect the normalization of the velocity and set Vp+ -^ vf- 

First we study (rf ^)stag at high temperature, i.e., T ^ 
TOp_. For the noninteracting case, the real part Xo be- 
comes logarithmic singular, i.e., x'oi'^^F + 'Zn,"', w, 2^) ~ 
{Attvf)~^ ln(A/a::) where x — Tii&x.{\vFq\\ |, |w|, T), while the 
imaginary part Xq is nonsingular The imaginary part of xo for 



12 



\uj\ ^ T is given by 

Xo(2A:i? + q\\,7T,u;,T) 



32vfT 



sech^(^). (4.25) 



By using the fact | Xq | S> | Xo I > the imaginary part of x is given 
by55 

x"{2kF + q\i,7r,u;,T) 
^Xstagl'?!!,'^,?') Xo(2A:f + g||,7r,w,r), (4.26) 

where x is an auxiliary function associated with the real part 
ofx: 

- ( T^-A dRex{2kp + q^^,TT,u;,T) 
Xstag(9lh^,7^) = 4^«F d\ii{A/x) ' ^ ^ ^ 

with X = niax(|cj|, \vFq\\\,T). The staggered part of the 
NMR relaxation rate can be obtained by inserting Eq. ( I4.26> 
into Eq. ( I4.18bt . Since Xoi'^kp + g||,7r,a;, T) has a sharp 
peak around g|| = with the exponential decay for large gy, 
we rewrite Xsta.giqh ^^ , T) ^^ XstaglT^Hn Eq. (I4.26> by ne- 
glecting the 511 and u; dependence in x^^ (g||,a;, T). Then 

the staggered part (Tf^) stag is given by 



(Tf^: 



stag 



SttvJ, 



vstag 



(T). 



(4.28) 



In order to examine Xsta g(^ )' we calculate the correlation 
function for n (x) [Eq. ( I4.4H which denotes the spin oper- 
ator with momentum transfer ^n w 2kp and q± — n. By 
applying the RG method to the spin-spin correlation function 

R{x,t) = (T.,n^(a;,r)n°:(0,0)), the RG equation for Xstag 
is obtained as 



^lnXstag(0 






1 



Q 



(4.29) 



at {Op^) = 
given by 



tt/ where / is integer Then the spin operators are 



n^ {x) ex cos 4>p+ sin 6c+ cos ^g.- , 
n^i (x) ex cos (/>p+ cos 9^+ cos (/ic- , 
nl (x) ex cos 0p_|_ cos 0CT+ sin 0^- 1 



(4.30a) 

(4.30b) 

(4.30c) 
where 0p+ is defined just after Eq. J4.8> . The RG equation for 

Xstag(^)isgivenby 



+ Mt + Ms 



^ inXstagl'J — 2 ^ 2 ''^ ^ 2 ^ 2 

(4.31) 
where the coupUng constants are calculated from Eq. ( I3.12t . 
Thus the NMR relaxation rate {T^^)sta.g in the region of 
niax{\mt\,ms) <C T ^ iTip_ is given by Eqs. J4.28I I and 
( 14.3 1> . However, in the present system, we find ttIs ~ rup-, 
and then this asymptotic behavior would not be realized. 

Next we calculate the antiferromagnetic contribution to the 
relaxation rate at low temperature by using the effective spin 
Hamiltonian ( 13. 7> . In this case, the staggered components of 
the spin operators ( I4.30> are rewritten as 



nl (x) w cos 0p+ (criAt2At3) ^4, 
n^ (x) « cos 4>p+ {^ll(72^J.S.) /^4, 

nl(a;) w cos ^p+ (^1/120-3) /14, 



(4.32a) 
(4.32b) 
(4.32c) 



where the Klein factors have been omitted. We follow the cal- 
culation performed in Refs. 58 67, and 68. If the spin-spin 
correlation function at zero temperature exhibits power-law 
behavior with exponent 77, i.e., (n_(a;)n_(0)) oc x^'', it can 
be shown that the temperature dependence of the NMR relax- 
ation rate is given by (Tj^^)stag oc T^~^, by using the con- 
is 



formal mapping techniquej^ If \mt\ 
realized, the exponent is 77 = ^-ftTp^ 



\ms\ <^ T <^ nip 
1 and then we obtain 



(Tf 



Jstag 



cxT 



K,+ /2 



(4.33) 



where I = ln(A/r) and Xsta, 
ature dependence of (Tj^^ 



(0) = 1. To obtain the temper- 
stag from Eq. (I4.28> . we first solve 



the RG equations (13. 1> for the coupling constants, and next 
substitute those into Eq. ( I4.29> . We note that, for the nonin- 
teracting case, the auxiliary function is given by Xstag(0 — 1 
and then (rj~^)stag = T/Snvp, which has the same tempera- 
ture dependence as the uniform part [Eq. ( I4.21H in the nonin- 
teracting limit. 



This behavior, however, would not be observed since m^ « 
■nip-. For \mt\ < T < |ms| in the SCd or CDWh-PDW 
phase, we can replace fi4 by its average value (/i4) 7^ in 
Eq. (I4.32> and then the exponent of the spin-spin correlation 



function becomes rj = ^K^ 



iT,-')s 



p+ 



cxT 



Thus we have 



-l/4+K,+ /2 



At temperature below the charge gap nip-, i.e., for I > In the limit of low temperature (T ^ \m.t\ 



(4.34) 



.), the 



Ip-, the field Op- 



can be replaced by its average value 



Without losing generality, we can assume that 9p- is locked 



relaxation rate exhibits thermally activated behavior. In the 
SCd phase [nit > 0), we obtain (see the Appendix) 
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(^1 ^)stag OC Y^ 



— cosh I ie—r] -\ i I 



2tt 



2 27rT' ' 



2TJ V vf 
2TrS(E{9i) - £{62) + E{03 



^) 



d9i (102 dOs 

'2^^^'8cosh[E{9i)/2T] cosh[E{e2)/2T] cosh[E{93)/2T] 

I 



coth^ 



(4.35) 
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FIG. 8: The temperature dependence of the NMR relaxation rate 

(rr^)stag for U/t = 3, Vn/t = V±/t = 0.7, t = t\\=t±^ 1, and 
(5 = 0.1. 



where E{9) = \mt\ cosh 6* is the rapidity representation of 
dispersion, B{x,y) is the beta function, and rj — Kp+/2. In 
the low temperature limit T -^ rrit, the staggered part of the 
NMR relaxation rate in the SCd state is given by (see the Ap- 
pendix) 



(rr')stag (X Ti+^^+Z^ exp(-2mt/T), 



(4.36) 



showing the activation behavior with a gap 2r7if . On the other 
hand, in the CDW+PDW state (m^ < 0), by considering two 
magnon process, the NMR relaxation rate would be given by 
(see the Appendix) 

(Tf i)stag « T^^+Z^ eM-\mt\/T) InT. (4.37) 



In Fig. |8] we show the staggered part of the NMR relax- 
ation rate, (Tj~^)stag at low temperature, which is calculated 
based on Eqs. J4.33> . ( I4.34> . and (I4.36> . The enhancement of 
{T^^)stag in the interval region of \mt\ < T < nis, which is 
similar to that of single chain, originates in the antiferromag- 
netic fluctuations. The rapid decrease of {T^^)sta.g at temper- 
ature of T < \mt\ shares the common feature with {T^^)uni 



andxs 



We note that the relative magnitude of {T^ )uiji and 
ag depends on those of their hyperfine couplings. 



C. Quantum critical behavior 

Finally we focus on the temperature dependence of Xs, 
(Tf^ )uni, and (Tf^ )st,ag, just above the QCP, i.e., m* = 0. 
Since the 3^+ term becomes marginally irrelevant in the ef- 
fective theory (13. 3> . we first examine the scaling of 5^+ which 
gives rise to logarithmic corrections for the physical quanti- 
ties. 

In the effective spin Hamiltonian ( 13. 3> . we consider the 
temperature region lower than the gap in the Majorana sin- 
glet sector. In this case, we can integrate out the ^^ degrees of 
freedom and then we can rewrite m° — i(7cr-(C+<^-) ^ "^t- 
Further we can set m^ = on the critical point. Then only 
the (7cr+ term remains where the RG equation is given by 



dl 



GMl) ^ -Gl+il), 



(4.38) 



and gcr+{l) ~ {2TTVp)Ga+{l)- By solving Eg. i43S\ . we have 



9<T+il) 



y<T+ 



l + {g^+/2TTVF){l-ls)' 



(4.39) 



where the initial value is given by ga-+{ls) — 9a+ = 
{27rvp) Gcr+{ls)- The quantity Ig corresponds to the scale of 
the gap in the Majorana singlet excitation, m^ w Ae~'=. 

From Eq. ( I4.38> . the character of the phase transition is de- 
termined by the sign of the initial value of Ga+ ?^-^^-^^ When 
the initial value is given by Ga+{ls) > 0, the coupling con- 
stant Gcr+{1) decreases to zero under renormalization, and 
then becomes marginally irrelevant. In this case, the effec- 
tive theory in the low-temperature limit leads to the noninter- 
acting massless Majorana fermion, and thus the system ex- 
hibits a quantum critical behavior. On the other hand, for 
the initial value being negative, i.e., Gcr+{ls) < 0, the cou- 
pling constant Ga+{1) becomes marginally relevant due to 
its divergence at It = 2tivf /\g(y+\- In this case, the effec- 
tive theory does not give a quantum critical behavior due to 
a mass gap rrit ~ Ae^^'^^^/I^^+I, even if the bare mass rn" 
reduces to zero. From the numerical calculation of Eqs. \'i.l\ 
and ( I3.12K we have confirmed that the coupling constant Ga+ 
ail — Is is positive within our choice of repulsive interactions, 
and that the present ladder system corresponds to the former 
case, i.e., the system exhibits the quantum critical behavior. 
If the sign of Ga+ih) could be changed by another type of 
interaction in the microscopic Hamiltoniani^iS the first-order 
transition would be obtained instead of the QCP within the 
present framework. 

In the following we examine the temperature dependence 
of Xs{T), (Tj"^)uni, and (Tj"^)stag, separately by using Eq. 
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FIG. 9: The temperature dependence of the NMR relaxation rate 
with f//i = 3, 1^1 /i = Vx/f = 0.7, t = t\\ = t^ ^ 1, and 
5 = 5c ~ 0.18. The power-law behavior is retained even in the 
limit of low temperature due to the vanishing of the Majorana triplet 
mass, i.e., rtit — 0. The dashed and dotted curves correspond to the 
contributions with q± — and n, respectively, while the solid curve 
denotes the total (Tj" )uni. 



j4.39> with / = ln(A/r) where A is a high-energy cutoff 
of the order of the bandwidth. We note that at low tempera- 
ture, i.e., for large /, the renormalized coupling constant shows 

g^+{l)^2TTVF/HA/T). 



1. The uniform spin susceptibility Xs{T) 



We can use the formula ( I4.13t for calculating the uniform 
spin susceptibility since rrit = at the QCP. By inserting Eq. 
( I4.39> into Eq. ( I4.13> , the temperature dependence of Xs (T) 
is obtained with the low-temperature asymptotics. 



XsiT) 



1 



2'kvf 



1 



ln(A/r) 



(4.40) 



2. The NMR relaxation rate: (Tj ^)uni and {T^ ^) 



stag 



The temperature dependence of (Tj^^)uni can be obtained 
by inserting Eq. (|^^J into Eq. ( 14. 2H . The second term of the 
RHS of Eq. ( 14. 2H can be discarded since this would show the 
exponential decay at temperature below rn^ as seen from the 
dotted curve of Fig.|5] For the low- temperature limit, we have 



{T^\ 



T 



Snvj, 



ln(A/T) 



(4.42) 



The overall temperature dependence of {T^ ^)uni at the QCP 
is shown in Fig.|9] For the single chain with S — ^, we have 



iT^'': 



ID 
uni 



T [xi^in 



Anvj, 



xliT) 



T 



invj, 



ln(A/T) 



(4.43) 



where the first equality is obtained in Ref.i55i. 

The logarithmic correction to the staggered part (T'j~^)stag 
can be obtained as follows. If one neglects the marginally 
irrelevant g^^, the staggered part (rj^^)stag is given by Eq. 
( I4.34> . In order to retain the renormalization effect of gcr+, 
we use Eq. ( I4.28> . In a way similar to the derivation of Eq. 
( I4.38> . the scaling equation of the auxiliary field Xsj^g (^) oi^ 
the QCP at temperature below nis is given by 



^lnxstag(0 = 4-^ 2 "^ 



+ ^G,+ (0, 



(4.44) 



where Gp+ is independent of I and the factor 3/4 in the RHS 
is determined from the scaling dimension of the spin operator 

Eq. ( I4.32> . i.e., (2 — 2dim[cos(/)p+] — 2dim[(Ti/i2/i3]) = 3/4. 



By solving Eq. (I4.44> with Kp == 1 — Gp+, we have 



AstagV / \ \ 



From Eqs. jA^W 
asymptotics of {T^ 



and 

1 



-5/4+Kp+/2 



(4.45) 
45}, we obtain the low-temperature 




J stag 



is obtained as 



Figure |6l shows the temperature dependence of Xs{T) where 
S = Sc corresponds to the QCP. Equation ( I4.40> is compared 
with the spin susceptibility for the S = ^ Heisenberg single 
chain given b}sS 



•^"<^'»2^ 



1 



1 



2 ln(A/T) 



(4.41) 



where v = TrJ/2 and J is the exchange interaction. The 
results ( I4.40> and ( I4.41> are consistent with the suscepti- 
bility in the 811(2)^ WZNW critical theory with marginally 
irrelevant operators, which gives the logarithmic correction 
{1 + fc/[2 ln(A/T)]} where k is the level of SU(2) algebra.™ 
Note that k = 2 for a two-leg ladder, while k — 1 for the 
S" = i single chain. 



m-': 



stag OC 



j.-l/4+K,W2 Ij, 



(4.46) 



This result is compared with the staggered component of T^ ^ 
for the single chain, which is given by 



(Tf') 



ID 
stag 



cxT 



K^ 



/In 



(4.47) 



We note that, in the insulating state {Kp -^ 0), Eq. ( I4.47> is 



y/\n{A/T) reproducing the result of 



reduced to (rfi)iDg 

the S = i Heisenberg spin chain,^' while Eq. ( 14.461 1 leads to 

(rf ^)stag -> T~'^/*y/\n{A/T) which is consistent with the 
result obtained in Ref. 58 except for the logarithmic correc- 
tion. 
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V. CONCLUSIONS AND DISCUSSION 

In the present paper, we have examined the ground-state 
phase diagram and the temperature dependence of the sus- 
ceptibility and the NMR relaxation rate for the extended 
two-leg Hubbard model away from half filling, by using the 
weak-coupling bosonization method. In the ground state, 
we have clarified the competition between the SCd state and 
the CDWh-PDW state and have shown the quantum criti- 
cal behavior close to the transition point where the SCd 
state changes into the CDWh-PDW state with increasing the 
nearest-neighbor repulsion and/or decreasing doping rate. At 
finite temperature, the magnetic response exhibits character- 
istic property coming from two modes of spin excitations. 
Especially on the quantum critical point, we found that the 
spin susceptibility shows paramagnetic temperature depen- 
dence with logarithmic corrections and the NMR relaxation 
rate exhibits anomalous power-law behavior. 

Here we discuss the commensurability effect due to the 
umklapp scattering which would play an important role close 
to half filling. At half filling, the umklapp scattering is given 
by^^3 

Humklapp ^lJ2Yl [^"y V'p,,.,^ ^laX2 ^-p,a,U '^-p.-.Ca 
P.o- Ci = ± 

(5.1) 

where 5^]^ = /,F_l +TO3,e^| andg^^ = {U + l,V±+m3^,V^i) 
with the numerical factors l± = ±1, TO3_+ = —1, and 
?7i3 _ = —2. In terms of bosonic fields, Eq. J5.1> is rewrit- 
ten as 



n 



1 



umklapp 



r[.9c 



cos2(/)p-|_ coa29p 



+ gc+,s+ cos20p+ cos2(l)cr+ 
+ gc+,s- cos2(/)p_(_ cos2(/)o._ 
+ .9c+,i^ cos2(/)p+ cos26cr-], (5.2) 



where the coupling constant are g^ 



-slw 



J3\\ ' 



-93± ' and 5^+^- 



^93± ' 9c+,s+ — 
: = +.g3j7- For 



the rung-singlet state at half filling, the renormalized coupling 
constants in Eq. ( 15. 2> are given by g* ^ < 0, g*_^ g_|_ < 0, 
g*^ g_ < 0, and g* -— ~ 0. Even in the presence of a 
few holes, the renormalized umklapp scattering would remain 
finite unless at extremely low energy scale. Then here we fix 
the amplitude of the umklapp scattering and discuss the effect 
of the finite doping. In terms of phase variable, the particle 
number operator is given by 



. c, , „ - 



2 



(5.3) 



From Eq. ( I5.3l l, the injection of a single electron or hole corre- 
sponds to the formation of the 7r/2 soli ton or antisoliton in the 



(f)p+ mode. In order to avoid the increase of energy, the 7r/2 
soliton in the 0p+ mode is always accompanied by the tt/2 
solitons in the 4'a+, dp-, and 0o-_ modes. This fact implies 
that the 7r/2 soliton in the <j)p^ mode involves the appearance 
of local spin at the same rung. On the other hand, the tt soliton 
in the 0p+ mode, which corresponds to iV = 2, is not accom- 
panied by solitons in the spin and other modes. Thus, if a tt/2 
soliton in the (/)p+ mode is created in the system, the free spin 
appears at the same rung in the bulk rung-singlet state and 
the spin-charge separation does not take place. This picture 
would connect with the strong-coupling one in the sense that 
holes can destroy spin singlets in the two-leg ladder systemsi^ 
Finally we compare the present results with the ex- 
perimental ones on the two-leg ladder compounds 
Sri4_a;Ca^Cu2404i, which have the characteristic fea- 
tures of the spin-gapped normal state and the superconducting 
state. For x = 12 and under a pressure of 3.5 GPa, the 
NMR measurements show two excitation modes above the 
SC state, where T^^ decreases rapidly at higher temperature 
and T-linear dependence is found at lower temperature.— 
This result resembles the present result of Figs. and |9] 
The decrease for T > iris comes from the formation of 
the spin gap in the spin-singlet excitations while the linear 
dependence for \mt\ < T < nis appears due to the gapless 
mode with the freedom of spin-triplet excitations. Such an 
interval region is enlarged close to the QCP since nit ~^ 
at the QCP, as seen from Fig. |3l From the NMR shifts, it 
is shown that the uniform magnetic susceptibility decreases 
slowly for 30 < T < 200 K and stays constant for T < 30 K. 
The former resembles Fig. |5] where the slow decrease below 
T/t < 0.3 is due to the band effect. The latter would 
correspond to Fig. |5] with S — 0.1, where Xs is almost 
independent of temperature for 0.07 < T/t < 0.01. Thus 
the present scenario of two spin excitations could be relevant 
to experiments although the present approach is based on 
weak-coupling theory. 
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APPENDIX A: STAGGERED PART OF THE NMR 
RELAXATION RATE IN THE LOW-TEMPERATURE LIMIT 



In this appendix, we derive the staggered part of the NMR 
relaxation rate (7\ ^)stag in the low-temperature limit [Eqs. 
J4.35> — ( r4.37H . based on the Majorana-fermion description of 
the effective theory. 

In the temperature region T ^ (l^^sl, m-p-), the staggered 
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component of the spin operator [Eq. (I4.32H is rewritten as 

nl oc cos 0P+ (0-1^2^3), (Ala) 

(Alb) 
(Ale) 



vL oc C0S(^p+ (/^i 0-2 /is), 
il ex cos^p+ (^1^20-3), 



where ^p_|_ = 0p+ — tt6x. In terms of these operators, the 
NMR relaxation rate is given by 



iT,-')s 



dt S{t), 



(A2) 



where S{t) — {nt{x = 0,t) n"(0, 0)) is the local correla- 
tion function at finite temperature. We will estimate this cor- 
relation function by using the effective Hamiltonian (13. 8> and 
(I4.24> . Since the charge and spin degrees of freedom are de- 
coupled, the correlation function can be rewritten as 



5(t) = 5p+(i) Rising (i), 



(A3) 



where 



Sp+{t) 

"Rising (''j 



(cos<^p+(0,i)cos</>p+(0,0)), (A4a) 
(/i(0,i)A*(0,0))2(a(0,i)a(0,0)).(A4b) 



One can easily find that all the correlation functions 
(n" (0,i)n" (0,0)) for a = x,y,z become identical, since 
the system has spin-rotational symmetry. 

The local correlation function for the charge fields is given 
byS 



Sp+{t) 



1 



-i(7r/2)?)sgn (t) 



■KTa/v 



sinh(7rr|<| 






2TiTa 



xB 



''2^ 



1-7? 



r;-l 



(A5) 



where r/ = Kp+/2 and B{x, y) is the beta function B{x, y) — 
T{x)r{y)/r{x + y). In the second equality, we have per- 
formed the Fourier transformation.*^ 

The correlation function for the Ising fields at finite tem- 
perature can be calculated following Ref. 72. The asymptotic 
behavior of the Ising correlation function depends on whether 
the system is in the ordered phase or in the disordered phase, 
i.e., depends on the sign of the mass nit. 

In the SCd phase (mt > 0), the Ising systems ai {i — 
1, 2, 3) are in the disordered phase; thus /i; has a nonzero ex- 
pectation value. In this case, the dominant contribution in the 
low-temperature limit isiSI, 



"-"Isingitjl 



T?lj>0 



^ r°° d0i de2 d03 

2-^ \ 97r 97r 97r 



^,^3 — ± 



2-K 2-K 2-K 



'ileE{6i)-eE{e2)+e3E{e3)]t 



X e 



(A6) 

where E{e) = \mt\cosh9 and f,{e) = [1 + g-^-^W/^]"!. 
By inserting Eqs. ( IA5> and ( IA6> into Eq. ( IA3> . and inserting 
it into Eq. ( IA2> . we obtain Eq. (I4.35> . 



In the CDWh-PDW phase (mt < 0), the Ising systems are 
now in the ordered phase; thus ctj has a nonzero expectation 
value. In this case the dominant contribution at low tempera- 
ture reads 



X M9,) /_,(02) e— [^(«0-B(92)]t_ (A7) 
Then the NMR relaxation rate for m^ < is given by 

duj 



(Ti l)stag oc Y^ 
r,-l 



cosh ( ie—T] + _ , 
27r \ 2 ' 2TJ 



IT LU \ 



X 



/27rT 

V VF 



^'2-^^2^' 



1-77 



dOi d92 2ttS{E{9i) - £{62) + lj) 

_oo ^^ 4cosh[£;(6li)/2T] cosh[^(6'2)/2r] ' 

(A8) 



The low-temperature asymptotic forms of (Tj^^)stag [Eqs. 
(I4.36> and ( I4.37H are obtained as follows. The u; integral in 
Eq. ( 14.351 1 or Eq. (JASjl is cut by the temperature ±T, since 
by performing the summation with respect to e {— ±) in Eq. 
"33 or Eq. (|A8} one obtains (for |tj| > T) 



yj cosh ( ie—T] + -rT= ) -B I -r — «e 



e=± 



2TJ 



sin(7r77) r(l — rj) 



(\^ 



\2ttT 



2^ 

r,-l 



1-r;) 



exp 1 -^ 



(A9) 



This can easily be verified by using the asymptotic form of the 
beta function: 

B{J--iS,\~ 7]) w r(l - 77) (-iS-)"^! (for S -^ ±00) 

(AlO) 
Thus we have for nit > 



(Tfi) 



stag 



OC 



T") / deid92de3 coth^ 



00 



X S{E{9i) ^ E{02) + E{03)) 

X g-[B(ei)+S(92)+-B(e3)]/2T 

oc Ti+"exp(-2mt/T), 



(All) 



which reproduces Eq. < I4.36^ . For mt < v/e have 

/oo 
deid02S{E{ei)-E{e2)) 
-oo 



X g-[s(ei)+B(e2)]/2T 
oc r''exp(-|mt|/T) InT, 

which reproduces Eq. ( I4.37> . 



(A12) 
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